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ABSTRACT 

This is the third of a series of papers of low X-ray luminosity galaxy clusters. In 
this work we present the weak lensing analysis of eight clusters, based on observations 
obtained with the Gemini Multi-Object Spectrograph in the g ', r' and i' passbands. 
For this purpose, we have developed a pipeline for the lensing analysis of ground-based 
images and we have performed tests applied to simulated data. 

We have determined the masses of seven galaxy clusters, six of them measured for the 
first time. For the four clusters with availably spectroscopic data, we find a general 
agreement between the velocity dispersions obtained via weak lensing assuming a 
Singular Isothermal Sphere profile, and those obtained from the redshift distribution 
of member galaxies. The correlation between our weak lensing mass determinations 
and the X-ray luminosities are suitably fitted by other observations of the M — Lx 
relation and models. 
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1 INTRODUCTION 


Clusters of galaxies are the most massive virialized struc¬ 
tures in the Universe. Hence, they are excellent laborato¬ 
ries to study the physics of baryonic and dark matter at 
large scales in bound objects llVoit 20051 : IPratt et~ahl 120091 ; 
lArnaud et al.ll201fll : iGiodini et a] .11201 ,'T . Numerical simula¬ 
tions show that massive clusters are formed from the mer¬ 
ging of smaller stru ctures in the hierarchical s tructure for¬ 
mation (see review, iKravtsov fc Borganil 1201 j ). Therefore, 
the study of low X-ray galaxy clusters could shed light on 
the assembly processes and environmental effects on their 
galaxy population, since these systems are likely to be e- 
volving by sub-structure interations and accretion. In these 
systems, velocity dispersions are lower than in massive clus- 
ter (< 800 km s^ 1 ), favoring the interactions and mergers 
between the galaxy members. Thus, morphological trans¬ 
formations are more frequent in these clusters. Also, low 
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mass clusters are more common than rich clusters due to the 
steepness of the cluster mass function. However, at the same 
time these systems are fainter and cooler, which makes them 
more difficult to detect and distinguish from background. 
Hence, these clusters have not been extensively studied com¬ 
pared to massive, luminous X-ray systems. 


The evolution of galaxy clusters has been probed to 
be determined by cosmological parameters. In particu¬ 
lar, the cluster mass function provides observational con¬ 
straints to cosmology, given its sensiti v ity on the cos¬ 


mological param e ters (e.g., Mandelbaum&_Selj|ds 
Rozo et al. 20091: VijihljMi etal.1 120091 : Allen et al 


2007 


2011 


Planck Collaboration et~ahl 20141 ). The main limitation in 


the use of this mass function is the practical deter¬ 
mination of the masses. Weak and strong gavitational 
lensing probe the projected mass distribution of clus¬ 
ters, with strong lensing confined to the central regions 
of clusters, whereas weak lensing can yield mass mea¬ 
surements for larger radii. Mass estimations from grav¬ 
itational lensing is affected by substructure, triaxiality, 
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large-scale structure and the possible presence o f mul 


Sereno 200?i; 

Corless & King 2009; Meneelretti et al. 

201C; 

Sereno et al. 

20ld: Sereno & Umetsu 2011 

: Giocoli et al. 

20121; Sereno & Zitrin 20121; Spinelli et al 

20121). 

How- 


ever, other methods such as the caustic technique em- 
ploying spectroscopic measurements of galaxies velocity 
mines fc Diaferic] 120061 ), might be expensive in telescope 
time. Besides, radial mass distribution of clusters could be 
determined using X-ray surface bright ness under t he as¬ 
sumption of hydrosta tic equilibrium (iLaRoaue et al.l|2006l : 
iDonahue et al.l 120141 '). Nevertheless, deviations from the 
equilibrium could highly affect the estimations. Therefore, 
gravitational lensing is an excellent and a fairly clean tech¬ 
nique for mass cluster determinations. 

Galaxy clusters and groups are expected to follow sim¬ 
ple relatio ns linking th e total mass with another physical 
quantities (lKaiseilll986i ). Given the difficulties of determin¬ 
ing the mass of these systems, the study of these relations are 
important since they are suitable to convert simple observ¬ 
ables into mass estimates. I 11 particular, the X-ray lumino¬ 
sity of groups and clusters can be considered a good tracer of 
halo m asses with approximately 20% scatter in the M-Lx re- 
lation jStanek_et i _al] 2006^ Maugha3||200^ Pratt_etjd. 2009|; 

stlViki 


IRozo etHh 20081 : Rvkoff et al. 20081 : Vikhlinin et alJ 20091) . 


The main advantage in its use is that X-ray luminosity can 
be accurately measured at high redshifts, requiring only 
previous cluster detection and redshift information. Weak 
lensing provides a suitable technique to study the M-Lx 


(Bardeau et al. 

20071: [Hoekstra 20071; Rvkoff et al.1 20081; 

Leauthaud et al 

2010: Okabe et al. 201(3). In this sense. 


three studies s panning from low X-ray luminosity clus¬ 
ters to groups (iRvkoff et al.l 120081 : iLeauthaud et all l20ld : 
iKettula et al.ll2014T) show a s ingle relation with a well de¬ 


fined slope ( Foex et al.ll20lj ). in agreement with those of 


massive clusters. 

This work is the third in a series of papers aimed to un¬ 
derstand the processes involved in the formation and evolu¬ 
tion of low X-ray luminosity galaxy clu sters at intermediate 
redsh ifts. The first paper of the series dNilo Castellon et al.l 
l20l4 hereafter Paper I) contains the main goals, sample 
selection, and details of observations and data reduction 
for both, photometry and sp ectroscopy. The second pa¬ 
per dNilo Castellon et al.1 12014. hereafter Paper II), presents 
photometric properties of seven low X-ray luminosity ob¬ 
served with Gemini telescopes. As the redshift increases, an 
increment of blue galaxies and a decline in the fraction of 
lenticulars is observed, while the early-type fraction remains 
almost constant. These results are in agreement with those 
for high mass clusters. At lower redshifts, the presence of 
a well-defined cluster red sequence extending by more than 
4 magnitudes showed that these intermediate mass clusters 
had reached a relaxed stage. 

In this oportunity we present the weak lensing analysis 
of eight galaxy clusters of the low X-ray luminosity sample. 
The paper is organizated as follow: In Sec. 12. 11 we describe 
the sample of clusters, and the acquisition and reduction of 
the images. In Sec. [3] we give the details of the weak lensing 
analysis for the mass determination. In Sec.|4]we present and 
discuss the estimated mass, and compared them with X-ray 
luminosity. Finally, in Sec.[5]we summarise the main results 


of this work. We adopt when necessary a standard cosmolog¬ 
ical model Ho = 70kms -1 Mpc -1 , =0.3, and LIa =0.7. 

2 GALAXY CLUSTERS, OBSERVATIONS AND 

DATA REDUCTION 

2.1 Sample description 

The studied sample of low X-ray luminosities was se- 
lected from the c atalogue of extended X-ray sources by 
iMulIis et al.l (120031 ). This catalogue is a revised version of the 
223 galaxy clusters serendipito usly detected in the R OSAT 
PSPC pointed observations by Ivikhlinin et ahl (Il998l ). Our 
galaxy cluster sample comprises a random selection of 19 
systems from the total sample of 140 galaxy clusters with 
X-ray luminosities in the [0.5-2.0] keV energy band (rest 
frame), close to the detection limit of the ROSAT PSPC sur¬ 
vey ranging from 10 42 to ~ 50 x 10 43 erg s -1 . The redshift 
range of our selection is 0.16 to 0.70 and a full description 
of the project and sample can be found in Paper I. 

The galaxy clusters subsample studied in this work is 
mainly based on the clusters optically analized in Paper II: 
7 galaxy clusters with X-ray luminosity ranging from 1.4 
to 26.1 xlO 43 ergs -1 in the [0.5-2.0] keV energy band, and 
redshifts between 0.185 to 0.7. We add to this sample with 
observed colours, the galaxy cluster [VMF98]102 located at 
z ~ 0.401, observed only in r' passband. In Table[T]we sum¬ 
marize the main characteristics of the clusters. The mean X- 
ray luminosity in [0.5-2.0] keV band is 13.4 x 10 43 erg s -1 
, an intermediate/low luminosity when compared to ~10 42 
erg s -1 for groups with extended X-ray emission or the larger 
values than 5 ~ 10 44 erg s -1 of rich clusters. Lx in [0.1 
2.4] keV band are used for further analysis and comparison 
with other workes (see Section[OJ 


2.2 Observations 

Photometric observations for the eight galaxy clusters were 
obtained with Gemini North (GN) and South (GS) tele¬ 
scopes, during the system verification process (SVP) and 
specific programs with Argentinian time allocation. Seven 
clusters w ere observed using the Gemini Multi-Object Spec¬ 
trograph (|Hook et all 120041 ) in the image mode, in the r' 
and g' or %' passbands with an array of three EEV CCDs 
of 2048 x 4608 pixels and only one (]VMF98]102) in the r' 
passband. Using a 2 x 2 binning, the pixel scale is 0.1454 arc- 
sec per pixel which corresponds to a FOV (Field of View ) 
of 5.5 x 5.5 arcmin 2 in the sky. 

All images were observed under excellent photometric 
conditions, with mean seeing values of 0.75, 0.66 and 0.74 
arcsec in the g' , r' and i! filters, respectively. Some observa¬ 
tions were made under exceptional weather conditions, such 
as those made to the galaxy cluster [VMF98]001, with a me¬ 
dian seeing of about 0.485 in the r' image. Further details 
about these observations are given in Paper II. Columns 6 
to 9 in Table[T]show a summary of the photometric obser¬ 
vations. 

All observations were p rocessed wit h the Gemini IRAF 
package vl.4 inside IRAlQ (lTodvlll993l ) . The images were 


1 IRAF is distributed by the National Optical Astronomy Obser- 
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Table 1 . Low X-ray luminosity Galaxy Cluster sample 


[VMF98] 

Id. 

a 

(J2000) 

8 

(J2000) 

Lx 

[0.5-2.0] keV 
(h/ 0 2 10 43 cgs) 

Lx 

[0.1-2.4] keV 
(h/ 0 2 10 43 cgs) 

z 

Program 

Id. 

g' 

r' 

i' 

001 

00 30 33.2 

+26 18 19 

26.1 

30.7 

0.500 

GN-2010B-Q-73 

- 

15x300 

15x150 

022 

02 06 23.4 

+15 11 16 

3.6 

3.8 

0.248 

GN-2003B-Q-10 

- 

4x300 

4x150 

093 

10 53 18.4 

+57 20 47 

1.4 

1.6 

0.340 

GN-2011A-Q-75 

- 

5x600 

4x150 

097 

11 17 26.1 

+07 43 35 

6.4 

7.7 

0.477 

GS-2003A-SV-206 

12x600 

7x900 

- 

102 

11 24 13.9 

-17 00 11 

8.1 

9.3 

0.407 

GS-2003A-SV-206 

- 

5x600 

- 

119 

12 21 24.5 

+49 18 13 

42.7 

53.6 

0.700 

GN-2011A-Q-75 

- 

7x190 

4x120 

124 

12 52 05.4 

-29 20 46 

3.4 

3.4 

0.188 

GS-2003A-SV-206 

5x300 

5x600 

- 

148 

13 42 49.1 

+40 28 11 

16.2 

21.4 

0.699 

GN-2011A-Q-75 

- 

7x190 

5x120 


Notes. Columns: (1), the cluster identificatio n; (2) and (3). t he equatorial coordinates of the X-ray centre; (4), the X-ray luminosity in 
the [0.5 - 2.0] keV energy band obtained from lYikhlinin et al.l d 1998Tl ; (5), shows the X-ray luminosity in t he [0.1 - 2 .4] keV ener gy band 
calculated using L 500 from the MCXC catalogue (Meta-C atalogue of X-ray Detected Clusters of Galaxies. iPiffaretti et ahll201 1''l: (6), 
the mean redshift for each cluster from lMullis et ahl j2003t ); (7), the Gemini Program identification; (8), (9) and (10), the number of 
exposures and individual exposure time in seconds for each passband. 


bias/overscan-subtracted, trimmed and flat-fielded. The fi¬ 
nal processed images were registered to a common pixel po¬ 
sition and then combined. 


3 WEAK LENSING ANALYSIS 

We developed a pipeline based on Python Language (ver¬ 
sion 2.7. Available at http://www.python.org) to make the 
lensing analysis. The pipeline computes the shear profile and 
fits a model to estimate the mass of a galaxy cluster, taking 
as input the observed image of the cluster. In the next sub¬ 
sections, we describe in detail the implemented weak lensing 
analysis pipeline and the results of the application on simu¬ 
lated data to test its performance. 

3.1 Object detection and classification 

The first step in the lensing analysis is the detection and 
classification of the sources in stars and galaxies. To perform 
the detectio n and photometry of th e sources we implement 
SExtractor (IBertin fe Arnoutal 19961 ). From SExtractor out¬ 
put, we use for the analysis the parameters: MAG-BEST, as 
the magnitude in each filter; MU-MAX, defined as the cen¬ 
tral surface brightness of the objects ( Umax ); FLUX_MAX 
as the peak flux above background; FWHM as the gaussian 
full width at half maximum; CLASS-STAR as the stellarity 
index and FLAG, which corresponds to the notes generated 
by SExtractor in the detection and measurement processes. 

SExtractor is run twice (in a two-pass mode): A first 
run is made to detect bright objects in order to estimate the 
seeing and the saturation level of each image, and a second 
run to do the final detection. The first run of SExtractor 
is made with a detection level of 5<r above the background. 
The seeing is estimated using the average FWHM of the 


vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 


point-like objects selected from the FWHM/MAG-BEST di¬ 
agram, since for these objects the FWHM is independent of 
the magnitude. Determining the seeing is important for the 
star-galaxy classification, given that SExtractor uses it to 
compute the stellarity index. The saturation level is esti¬ 
mated as 0.8 times the maximum value of the FLUX_MAX 
parameter. These parameters, seeing and saturation level , 
are taken into account in the SExtractor configuration file 
for the second run, with a lower threshold detection limit of 
1.5cr. A second run is made in dual mode, detecting objects 
on the r' image, while astrometric and photometric parame¬ 
ters are measured on all individual images. 

Sources are classified according to their position in the mag¬ 
nitude/central flux diagram, the FW HM respect to the see - 
ing and the stellarity index, following iBardeau et ali (|2005l ). 
in stars, galaxies and false detections. In Figure[l]we show, 
as an example, Umax as a function of the r' magnitude 
(upper panel) and the r' magnitude against the FWHM 
(lower panel), for all objects dected by SExtractor in the 
cluster [VMF098]102. Objects that are more sharply peaked 
than the Point Spread Function (from now on PSF), thus 
with FWHM < seeing - 0.5 pixel, and with FLAG param¬ 
eter > 4, are considered as false detections. As the light 
distribution of a point source scales with magnitude, ob¬ 
jects on the line magnitude/central flux, ± 0.4 magnitudes, 
FWHM < seeing + 1 pixel and CLASS-STAR > 0.8 are con¬ 
sidered as stars. The rest of the objects are considered as 
galaxies. 

The first step in the pipeline ends generating two catalogues, 
one for the objects classified as stars and another for the 
galaxies. 


3.2 Shape measurements 

Measurements of galaxy shape are central in this analysis, 
given that galaxy ellipticities are used for the shear estima¬ 
tions and therefore to estimate cluster masses. It is impor¬ 
tant to take into account the roundness effects of the at¬ 
mosphere as well as the distortions caused by the telescope 
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(8 1 ,8 2 ) in the image plane, using a matrix transformation: 


S/3 1 

60 2 


1 — k — 71 
—72 


—72 

1 — K + 71 


SB 1 

SB 2 


where 71 and 72 are the components of the complex shear 
7 = 71 + * 72 - This can also be expressed as: 


5/3 1 

S/3 2 


**(!-«) 


1 - 51 -52 

-52 1 + gi 

where gn and <72 are the components of the reduced shear: 


SB 1 

SB 2 


0 = 


7 


1 — K 


(i) 


which is a nonlinear function of the two lensing functions: 
the complex shear, 7 , and the convergence, k, related to 
the projected mass density. If lensing is weak, the image of 
a circular source with ratio r, appears elliptical, with axis 
given by 


1 — « - |7| ’ 

Defining the ellipticity as 
a — b 


b = 


i-« + H 


__ l7l 

a + b 1 — k 


M 


where g becomes the normal shear, 7 , since k 1 , which 
generally holds in the weak lensing regime for clusters, and 
will be assumed henceforth here. 

If the source has an intrinsic ellipticity e s , the observed el¬ 
lipticity in the weak lensing limit will be: 


Figure 1. Classification of objects detected in the r' image of 
the galaxy cluster [VMF98J102. Here stars are represented by tri¬ 
angles, galaxies by points, and artifacts by cruxes. Upper pannel 
shows fiMAX / t' plane, where stars are situated in the region 
marked by the solid line ± 0.4 magnitudes, and in the lower pan¬ 
nel we show r'/FWHM plane. 

optics, all together included in the PSF, which is convolved 
with the galaxy intensity distribution. 

For the shape measurements we use IM2SHAPE 
(IBridle et al.l 12002 1. This code computes the shape param¬ 
eters modeling the object as a mixture of Gaussians, con¬ 
volved with a PSF which is also a sum of Gaussians. For 
simplicity both, the PSF and the object, are modeled with 
a single elliptical Gaussian profile. 

The PSF field across the image is estimated from the 
shape of the stars, since they are intrinsically point-like 
objects. We only used objects with a measured ellipticity 
smaller than 0.2 to remove most of the remaining false de¬ 
tections and faint galaxies present in the catalogue. Looking 
at the 5 nearest stars at each position, we have also removed 
those that differ by more than 2 cr from the local average 
shape. Then, we linearly interpolate the local PSF at each 
object position by averaging the shapes of the five closest 
stars. After PSF determination, we use again IM2SHAPE 
to measure the galaxy shapes, and the result is a catalogue 
of the galaxies with its intrinsic shape parameters. 

3.3 Shear radial profiles 

Gravitational lensing maps the unlensed image in the source 
plane, specified by coordinates (/I 1 ,/? 2 ), to the lensed image 


e = e s +7 

Assuming that unlensed galaxies are randomly oriented on 
the sky plane ((e s ) = 0 ) and averaging over sufficiently 
many sources: 

<e> = ( 7 ) (2) 

Hence, in the weak lensing approximation, we get an unbi¬ 
ased estimator of the reduced shear by averaging the shape 
of background galaxies in concentric annuli around the clus¬ 
ter centre. Spherical symmetry also implies that the average 
in annular bins of the tangential component ellipticity of the 
lensed galaxies, defined as the E-mode, traces the reduced 
shear. On the other hand, the average in annular bins of the 
component tilted at 7 r /4 relative to the tangential compo¬ 
nent, the B-mode, shoul d be exactl y zero for the case of per¬ 
fect symmetry (e.g. iBartelmann fe Schneidej|200ll . Sec. 4). 

Because of the random orientation of the galaxies in 
the source plane, the error in the observed galaxy clliptici- 
ties and thus, on the estimated shear, will depend on the 
number o f g alax i es ave raged together to measure the shear 
dSchneider et al.l l2000l ). Thus, the errors in the measured 
shear can be estimated as: 



where a e is the dispersion of the intrinsic ellipticity distri¬ 
bution (cr e « 0.3) and N is the number of objects in the 
annular bin. 

We have adopted the brightest cluster galaxy in r' filter 
as the cluster centre, a c riterium commonly used for lensing 
masses determination s llOkabe et al.l l20ld : iHoekstra et al.l 
2011; Focx et al. 2012). Shear profiles were computed using 
nonoverlapping logarithmic annuli, in order to have similar 
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signal-to-noise ratio (S/N) in each annuli. We have tested 
different annuli sizes but the final mass results have not 
showed a strong dependence on this parameter. We have 
fixed the size for the one we obtained lowest errors for the 
SIS and NFW profile fits. The profiles were fitted from the 
inner part were the signal becomes significantly positive, 
to reduce the impact of miscentering, up to the bin with 
highest number of galaxies (~ 3 arcmin for most of the clus¬ 
ters, which roughly corresponds to 0.8-1.4 Mpc). Our profiles 
were mainly limited by the FOV of the images. With these 
limits, 4-6 points were available in the shear profiles. 

3.4 Background Galaxies selection and redshift 
distribution 

To perform the shear estimation, background galaxies were 
selected as those with r' magnitudes between mp and 
rrimax + 0.5. mp is defined as the faintest magnitude where 
the probability that the galaxy is behind the cluster is higher 
than 0.7 and m maa; corresponds to the peak of the magni¬ 
tude distribution of galaxies in the r' passband. Keeping 
galaxies brighter than mmax + 0.5 ensures that we are not 
taking into account too faint galaxies with higher uncertain¬ 
ties in the shape measurements. We have also restricted the 
objects to those with good S/N and with a good pixel sam¬ 
pling by using only the galaxies with <r e < 0.2 (cr e is defined 
as the quadratic sum of the errors cr e i and cr e 2 given by 
IM2SHAPE) and with FWHM > 5 pixels. 

Once we obtain a catalogue for the background gala¬ 
xies, we average the components of the ellipticities (E-mode 
and B-mode) in nonoverlapping annuli. The average E-mode 
components corresponds to the shear value which depends 
on the geometrical factor fi = Dls/Ds, where Dls is the 
angular diameter distance from the lens to the background 
source galaxy, and Ds is the distance from the observer to 
the background galaxy. A galaxy at the same radial distance 
from the centre of the cluster but at a different background 
redshift is sheared differently. This variation is taken into 
account once we fit the profiles by (/3). 

To estimate mp and (/3) we used the catal ogue of 
photometric redshifts computed by ICoupon et al.1 (|2009h . 
based on the public release Deep Field 1 of the Canada- 
France-Hawaii Telescope Legacy Survey, which is complete 
up to m r = 26.. We compute the fraction of galaxies with 
2 > ^cluster in magnitude bins of 0.25 magnitudes for the 
r' filter, and then we chose mp as the lowest magnitude for 
which the fraction of galaxies was greater than 0.7. Then we 
applied the photometric selection criteria to the catalogue 
(mp < m r < m,nax + 0.5) and we computed fi for the whole 
distribution of galaxies. To take into account the contamina¬ 
tion by foreground galaxies given our selection criteria, we 
set j3{z p hot < Zciust.er) = 0 which outbalances the dilution 
of the shear signal by these unlensed galaxies. Deep Field 1 
covers a sky region of 1 degree 2 , thus to estimate the cosmic 
variance, we divide the field in 25 non-overlapping areas of 
~ 140 arcmin 2 and we compute mp and (/3) at z c i us ter = 0.5 
for each area. The uncertainties due to the cosmic variance 
were estimated as the dispersion of the values obtained for 
each area, obtaining ~ 0.3 for mp and ~ 0.01 for (/3). Given 
that the errors in (/?) are lower than the 3%, which represents 
an error of the ~ 5% in the mass, we did not consider these 
uncertanties in the estimation of the masses errors since the 



Figure 2. Fraction of galaxies with z > 0.5 (JV(z > 0 .5)/Ntot), 
for different magnitudes in filter r' and colours r'- i \ com puted 
using photometric redshifts given by ICoupon et ahl li2009l l, used 
to compute the weight for the shear estimation. The vertical line 
idicates the m max position (see text for its definition). 

uncertainty due to the intrinsic shape of field galaxies is 
much bigger. 

In order to take into account the contamination of fore¬ 
ground galaxies in the catalogue, we weighted the estimated 
shear, ( 7 ), with the probability that the gala xy was b ehind 
the cluster. We compute this probability usina lCoupon et al.l 
catalogue, from the fraction of galaxies with z > z c i us ter for 
each bin in magnitude, r', and colour (g’ - r’ and r’ - i’), 
see Figure[2] Hence, given the magnitude and the colour of 
each galaxy, we assigned to it a weigh, w, as the fraction of 
galaxies with z > z c i us ter in that bin. For [VMF98]102 we 
have only one image in the filter r', therefore for weighing 
the shear profile we take into account the probability that 
each galaxy was behind the cluster given the magnitude of 
that galaxy (we did not take into account the colours for 
computing this probability, as in the other clusters). 

3.5 Fitting the profiles 

We finally estimate the M 200 mass, defined as 
M 200 =M ( < R 200 ) = 200 p cr it(z)- 7 rr 2 oo, where R 200 
is the radius that encloses a mean density equal to 200 
times the critical density ( p cr it = 3H 2 (z)/8nG\ H(z) is the 
redshift dependent Hubble parameter and G is the gravita¬ 
tional constant). In order to do that we fit the shear data 
with t he singular isotherm al sphere (SIS) and the NFW 
profile ([Navarro et al.l 119971 ) using \ 2 minimization. These 
density profiles are the standard parametric models used 
in lensing analysis to characterize the lenses. Following, we 
explain briefly the lensing formulae for these two profiles: 

3.5.1 SIS profile 

The SIS mass model is the simplest one for describing a 
relaxed massive sphere with a constant and isotropic ve¬ 
locity dispersion. This is mainly described by the density 
distribution: 


This model corresponds to a distribution of self-gravitating 
particles where the velocity distribution at all radii is a 
Maxwellian with one dimensional velocity dispersion, ov. 





















6 Gonzalez et al. 


Semi-major axis map before PSF correction 



Semi-major axis map after PSF correction 



Figure 3. PSF treatment applied to stars of one of the images of the DES simulation: Mayor semiaxis (a cos 0 , a sin 6) before (left) 
and after (right) the PSF deconvolution in the CCD. Notice that the semiaxis are more randomly distributed and the scale (given by 
the first thicker segment in the upper-left corner and which corresponds to 3 pix) is much more smaller after the taking into account the 
PSF. 
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Figure 4. Shear profiles obtained for the Low (left) and High (right) Noise PSF Applied File, from the DES cluster simulation. The 
dashed curve shows the SIS profile for the input value of cry and the solid one the fitted profile. E and B modes are represented by full 
circles and crosses, respectively. 


From this equation, we can get the critical Einstein radius 
for the source sample as: 


Oe 


47TCT^ 1 

C lel (P) 


(4) 


where c ve i is the speed of light, in terms of wich one obtains: 


K0 = 79 = 


Oe 

20 


(5) 


where 0 is the distance to the cluster centre. Hence, fitting 
the shear for different radius, we can estimate the Einstein 
radius, and fro m that, we can obtain an estimation of the 
mass M 200 as Leonard fe Kindl2010l ~): 


M 2 00 = 


2 aj 

\/50 GH{z) 


(6) 


3.5.2 NFW profile 

The NFW profile is derived from fitting the density pro¬ 
file of numerical simulations of cold dark matter halos 


(jNavarro et al.Hl997h . This profile depends on two param¬ 
eters, the virial radius, R 200 , and a dimensionless concen¬ 
tration parameter, c: 

/ \ _ pc$c 

P T (r/r s )(l + r/r s ) 2 

where r s is the scale radius, r s =R, 2 oo/c and 5 C is the cha¬ 
racteristic overdensity of the halo, 

_ 200 __ 

3 ln(l + c) — c/(l + c) 

We used th e lensing formulae fo r the s pherical NFW density 
profile from I Wright fe Brainerdl (|2000T ). If we fit the shear for 
different radius we can have an estimation of the parame¬ 
ters c and R, 2 oo- Once we obtain R 200 we can compute the 
M 200 mass. Nevertheless, there is a well-known degeneracy 
between the parameters R 200 and c when fitting the shear 
profile in the weak lensing regime. This is due to the lack of 
information on the mass distribution near the cluster cen¬ 
tre and only a combination of strong and weak lensing can 
raise it and provide useful constraints on the concentration 
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parameter. Since we do not have strong lensing modeling for 
the clusters in the sample, we decided to fix the concentra¬ 
tion parameter, C 20 0 = 4, accor ding to the predicted con¬ 
centrations given by I Duffy et al.l (l201lh for a relaxed cluster 
with M= 1 x 10 14 Mq/iy O 1 placed at 2 ~ 0.4. Thus, we fit the 
mass profile with only one free parameter, R, 2 oo- 


[6] with the reduced x* f° r each fit. We include both fits, 
SIS (solid line) and NFW (dashed line) models. Points and 
crosses represent the E and B modes averaged in annular 
bins, respectively. All profiles are well fitted by both mo¬ 
dels. In the next subsections we discuss our results and we 
study the relation between the mass derived and the cluster 
X-ray luminosities. 


3.6 Testing the pipeline with simulated data 

To check the performance of our weak lensing analysis 
pipeline, we tested it on the D ES clu ster simulation images 
publically available dGill et al.l 120091 '). This simulation con¬ 
sists of a sets of images, with different grades of difficulty, of 
sheared galaxies due to the presence of a SIS profile with a 
velocity dispersion of 1250kms _1 . This is a suitable test for 
our pipeline given that the idea is to apply it to real clusters 
of galaxies. We applied our pipeline to three of the available 
image files, High Noise File, High Noise PSF Applied File 
and Low Noise PSF Applied File. 

For the PSF Applied files, we checked that our 
IM2SHAPE implementation can recover point-like objects 
by applying the PSF correction to each star. Figure[3] shows 
the results of the shape parameters measurements for these 
stars, with and wihtout taking into account the PSF in the 
shape measurement: the size distribution is dominated by 
point sources, and the orientation is more uniformly dis¬ 
tributed after the PSF correction. 

The images contain only the sheared galaxies, hence 
all the galaxies detected were considered as background 
galaxies at 2 = 0.8, which is the average redshift of the 
galaxies. We cut the catalogue discarding the galaxies with 
FWHM < 5 and with a e > 0.2. Shear profiles are shown in 
Figurc[4] For the most complex image that we treated (high 
noise image of sheared galaxies convolved with a PSF), we 
obtained a deviation parameter of 1.3, defined as the num¬ 
ber of a that the result is away from the input value of 

, , . result —input 

cry = 1250 km/sec, i.e. a = - , where the e- 

error 

rror was estimated according to the root mean square error 
of the Einstein radius. Given these results, we conclude that 
our weak lensing pipeline is able to reproduce the input shear 
signal, thus it could be applied to real observations to ex¬ 
tract the lensing signal and to estimate the masses of cluster 
of galaxies. 


4 RESULTS 


4.1 Properties of individual clusters 

4.1.1 [VMFJ001 

For the galaxy cluster [VMF98]001 we obtained a shear sig¬ 
nal consistent with a velocity dispersion of ^SOOkrns^ 1 . 
There is a big offset between the position BCG, adopted 
as the centre for the lensing analysis, and the X-ray lumi¬ 
nosity peak from ROSAT (~110kpc), not observed in the 
X-ray countours obtained with XMM-Newton (see Figure 6, 
from Paper I). Thus, given the lower resolution of ROSAT 
observations, the X-ray peak might be poorly determinated 
leading to unrealistic offsets. Further evidence of this fact 
is the absence of the shear profile signal centred at the X- 
ray position. We argue that the centre of the gravitational 
potential should be close to the BCG position. 

4.1.2 [VMFJ022 

The galaxy cluster [VMF98]022 shows an elongated distri¬ 
bution of galaxies in the NE-SW direction. The cluster is 
dominated by a bright elliptical galaxy, which presents a 
shift of ~12" in the south-west direction with respect to the 
X-ray peak emission (for further details about the cluster 
morphology, see Section 3.4 in Paper II). We compute the 
shear profile taking this bright elliptical as the centre of the 
cluster. This system presents a shear profile signal consistent 
with a velocity dispersion of 540 km s, in good agreement 
with the velocity dispersion fitted from the redshift distri¬ 
bution (see Sec. 4 in PaperI). 

4.1.3 [VMFJ093 

For the cluster [VMF98]093, in spite of the low density of 
background galaxies, we obtain a significant signal consis¬ 
tent with a velocity dispersion of 750 kms -1 . As evidence of 
the relaxed state of this cluster, we observe a dominant pop¬ 
ulation of red galaxies as well as concentric X-ray countours 
centred in the BCG (Figure 6 in Paper I). 


From our weak lensing analysis, we estimated the mass of 
seven clusters in the sample. Due to its low signal to noise, 
for cluster [VMF98]148 it was not possible to derive a re¬ 
liable mass estimate from our lensing measurements. The 
results of the analysis are shown in Table[2] Errors in cry, 
R 200 and the masses were computed according to the \ 2 dis¬ 
persion. Errors in are higher than Mfojf, given the 

big uncertainties in the R 200 parameter . Nevertheless, both 
estimations are consistent being the NFW masses systema¬ 
tically larger by a ~ 20% (( / Mf 0 J 0 s ) = 1.21 ± 0.13, 

where the uncertainty corresponds to the scatter around the 
mean), i n ex cellent agreement with the result presented by 
lOkabe et al.l (l201di l for the virial masses. Shear profiles ob¬ 
tained for the galaxy clusters are shown in Figures[5] and 


4.1.4 [VMFJ097 

The g ala xy clu s ter [V MF981097 was previously analysed by 
ICarrasco et al.l (j2007i l. using the same set of images. They 
obtained a large discrepancy between mass estimates, where 
the X-ray mass exceeds by more than a factor three the 
weak lens derived estimate. Moreover, they found a large 
degree of substructure, as also seen in the redshift distribu¬ 
tion presented in Paper I (Figure 11). However, substruc¬ 
ture cannot explain the defect in the weak lensing mass, 
given that substructure in the surroundings would tend to 
dilute the tangential shear leading to mass under-e stimation 
llMeneghetti et al]|20ld : iGiocoli et al.ll201 2|, |2014|) . We im- 
prove the profile (see Figure 9 from Carrasco et al.l 120071 . p. 
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r [Mpc h-, 0 1 ] 
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3 . \ [VMF98]001 |_ 

\ T - - Fit NFW profile R 200 = 1.3 ± 0.2 Mpc , x r 2 ed =1.2 

„ \ — Fit SIS profile a = 780 ± 100 km/s, x r 2 ed =0.8 
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Figure 5. Shear radial profiles as a function of cluster-centric projected distance (in arcsec and Mpc) obtained for the r' images of sample 
of clusters. The solid and the dashed lines represent the best fit of SIS and NFW profiles, respectively, with the fitted parameters given 
in the box. The points and crossings show the E and B modes profiles averaged in annular bins, respectively. Error bars are computed 
according to Eauationl3l 
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Table 2. Main results of the weak lensing analysis 


[ VMF 98 ] 

Id. 

a 

( J 2000 ) 

<5 

( J 2000 ) 

Pback 

mp 

TLlmax 

</ 3 > 

spec 

°v 

SIS 

a v 

M200 

NFW 

R200 M200 

001 

00 30 34.0 

+26 18 10 

56 

23.0 

26.1 

0.41 

- 

780 ± 100 

3.4 ± 1.3 

1 q + 0.2 
±0 - 0.2 

4 0"*” 2,2 
^• u - 2.0 

022 

02 06 21.2 

+15 11 01 

18 

20.7 

25.1 

0.61 

508 

570 ± 100 

1.5 ± 0.8 

-1 1 + 0.2 
1,1 - 0.2 

2 -i t\i 

093 

10 53 18.9 

+57 20 45 

8 

22.3 

24.0 

0.48 

- 

750 ± 140 

3.4 ± 1.9 

1 4 +o- 4 
-0.4 


097 

11 17 26.1 

+07 43 35 

40 

23.0 

26.0 

0.43 

775 

720 ± 100 

2.7 ± 1.1 

1 1 +0.3 
1 - ± -0.2 

2 . 8 +? 

102 

11 24 05.8 

-17 00 50 

40 

22.7 

25.9 

0.49 

675 

650 ± 120 

2.1 ± 1.2 

1 9 +0.3 
-0.2 

9 7 + 1-9 

119 

12 21 29.3 

+49 18 40 

13 

24.5 

25.4 

0.29 

- 

1000 ± 160 

6.3 ± 3.1 

1-4 ! g ;2 

7 q+3.8 
'•^-3.4 

124 

12 52 04.1 

-29 20 29 

33 

19.5 

25.7 

0.71 

700 

430 ± 60 

0.7 ± 0.3 

0.8 1 °;® 

o.sll;? 

148 

- 

- 

26 

24.5 

25.9 

0.29 

- 

- 

- 

- 


Notes. Columns: (1) shows the cluster identification; (2) and (3), the coordinates of the centre adopted for the lensing analysis; (4), 
the density of background galaxies (galaxies arcmin~ 2 ); (5) and (6), the brightest and faintest magnitude limits considered for the 
galaxy background selection (see Section ^. 411 : (7), the geometrical factor; (8), the line-of-sight spectroscopic velocity dispersion from 
Paperl; (9) and (10) the results from the SIS profile fit, the velocity dispersion and M 200 (see Eauationsl4l and l6)l : (11) and (12), the 
results from the NFW profile fit, R 200 and M 200 . The velocity dispersion, M 200 and R 200 are in units of kms 1 ,10 1 ' 1 MqK^q and 
Mpc hjQ , respectively. 


r [Mpc h 7 I ] 
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Figure 6. Shear radial profiles as a function of cluster-centric projected distance (in arcsec and Mpc) obtained for the r' images of sample 
of clusters. The solid and the dashed lines represent the best fit of SIS and NFW profiles, respectively, with the fitted parameters given 
in the box. The points and crossings show the E and B modes profiles averaged in annular bins, respectively. Error bars are computed 
according to Equation[3] 


11), adding a new constraint for the E-mode and obtain a 
profile consistent with zero for the B-mode. Nevertheless, 
our weak lensing mass estimat e is consistent with that ob¬ 
tained by lCarrasco et aD (I 2 OO 7 I ). corresponding to a velocity 
dispersion of ~ 700 km s _1 and also, with the velocity disper¬ 
sion from the redshift distribution of 775kms _1 (see Sec. 4 
in Paperl). 


4.1.5 [VMFJ102 

The results from cluster [VMF98]102 give a velocity disper¬ 
sion of 640kms _1 , in good agreement with the spectro¬ 
scopic value obtained in Paper I. In this case, the profile 
was built adopting a centre between the X-ray peak and 
the second brightest galaxy member. This was selected af¬ 
ter trying to fit the profile taken the centre as the second 
brightest galaxy and then, as the X-ray peak, without ge¬ 
tting enough signal-noise ratio to fit the profile. The sec- 
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log 10 (L x E(z) 1 [ergs 1 h 70 2 ]) 


Figure 7. Weak tensing m asses versus X-ray luminosities for the sam ple of clusters (d iamonds), combined with the stacked measurement 
bv TLeauthaud et ah feoich (ope n triangles), EXCPRES clusters by [Foex et al.l d2012h (open circle s) and low-mass f rom the CFHTLS 
(open squares) bv|Kettula et al] d2014h ■ Dashed, pointed, and solid lines represent the fit obained bv lLeauthaud et alJ (l201()l ). lFoex et al.l 
d2012h and iKettula et alTil^oTdl) . respectively. 


ond brightest galaxy in this case is close to the other bright 
cluster members and, unlike the brightest galaxy, it is an 
elliptical galaxy, so it is a more adequate guess for the clus¬ 
ter centre in this case. This cluster presents irregular X-ray 
contours (Figure 6 in Paper I) and, based in spectroscopic 
information, we found a non related group of galaxies in the 
line of sight (Figure 11, in Paper I). Also, there is a big off¬ 
set between the X-ray peak and the centre adopted for the 
lensing analysis (~220kpc), however we could not confirm 
this offset with higher resolution observations. 


agrees with the BCG position. Besides, there is no evidence 
of another group in redshift space (see Figure 11 from Pa¬ 
per I) and we observed a dominant red galaxy population 
(see Paper II), which indicates the relaxed state of this sys¬ 
tem. This cluster presents a low shear signal consistent with 
430 kms -1 . There is a large difference between the velocity 
dispersion obtained by the lensing analysis and that derived 
from the redshift distribution (700kms -1 , Sec. 4 in Paper 
I). We notice, however, the high uncertainty in this value 
given the small number of available redshifts. 


4.1.6 [VMFJ119 

The cluster [VMF98]119 is one of the highest redshift clus¬ 
ters (z ~ 0.7) in our sample. Even with a very low density of 
background galaxies, it shows a significant shear signal ac¬ 
cording with a velocity dispersion of 1000 kms -1 . The centre 
was placed at the brightest galaxy member, ~ 1’ from the 
ROSAT X-ray peak. Using X-ray observations from CHAN¬ 
DRA Data Archive, we built the X-ray contours and the 
peak is displaced ~0.9’ from the ROSAT centre, but still 
~0.4’ (~170kpc) displaced from the BCG. Also, the B- 
modes do not follow a null flat profile, which could be sug¬ 
gesting a large deviation from the spherical symmetry. This 
can also be seen in the distribution of member galaxies (Fig¬ 
ure 12 in Paper II). 


4.1.7 [VMFJ124 

Finally, for the cluster [VMF98] 124, the centre from the X- 
ray data using XMM-Newton contours (Figure 6 in Paper I) 


4.2 M — Lx relation 


We have also investigated the relation between the es¬ 
timated mass and the X-ray luminosity, which is a di¬ 
agnostic of the halo baryon fraction_and the entropy 
structure of the intracluster gas (iRvkoff et al.ll2008l) . The 
L x — M relation has been extensively studied , main ly 
at low redshift,s (z < 0.11 nsiriu Y-rav data, (Ma.rkevitch 


1998|; Arnaud et al 

200^; Reiprich & Bohrineerl 

20021; 

PoDesso et al. 200(3: 

Morandi et al. 20071: Pratt et al. 

2009; 

Vikhlinin et al. 2009 

). The main conclusion was that the 


relation follows a power-law, but with a slope and ampli¬ 
tude that differ from the self-similar prediction of M oc 
L^- 4 . Instead, they found a flatter slope, a = 0.56 — 0.63. 
Physical mechanism ruling the baryonic content of clusters, 
could strongly affect the X-ray luminosity, and so on the 
L x — M relation, causing deviations from a simple gravita¬ 
tional model. Simulations combining the gravitational evo¬ 
lution of dark matter structures together with the hydrody- 
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nami cal behaviour feorgani et al.ll2004l ; lKavll20()3 : iBorganil 
l2008h favor a lower slope value. 

Figure[7] shows the M — Lx relation for the galaxy 
clusters studied in this paper with masses estimated from 
the weak lensing analysis, together with those derived by 
other studies, the M — Lx rel ation based on 12 lo w mass 
clusters from the CFHTLS by iKettula et al] (12014) ; 11 X- 
ray bright clusters select ed and 206 stacked gala xy groups 
in the COSMOS held bv ILeauthaud et al.l d201Ch . and the 
Lx ~ M re l ation obtained from the EXCPRES sample by 
iFoex et ~ahl d2012h . In principle, the slopes from M — Lx 
(f3 ) and Lx — M (a) could be easily compared (a = 1 /ft) 
assu ming that the halo m ass function is locally a power- 
law (ILeauthaud et alj|2010f) . For comparison with other au¬ 
thors estimates, we used the N FW m asses showed in Ta- 
blc[?l Given that IKettula et ali (I20l4 derived core-excised 
luminosities, they are systematically lower than the rest of 
the plotted luminosities for a given mass. We notice that 
our mass determ i nation s are in very good agreement with 
ILeauthaud et all d2Qloli fit. The largest deviation from this 
fit corresponds to the two lowest X-ray luminosity clusters. 
Besides, [VMF98]093 contains a very low density of back¬ 
ground galaxies which affects the precision of the shear es¬ 
timates, and in the field of [VMF98] 124 there is a star with 
X-ray emition, which could bias high the quoted X-ray lu¬ 
minosity of the cluster. 


5 SUMMARY AND CONCLUSIONS 

In this work we presented the weak lensing analysis of 
eight low X-ray luminosity galaxy cluster. We described 
the pipeline for determining weak lensing masses of clusters 
using ground-based images. The analysis consisted in: the 
detection and classification of the sources, the shape mea¬ 
surements on the r' images taking into account the PSF, 
the galaxy background selection, the computation of shear 
profiles weighing the cllipticities according to the r' magni¬ 
tude and the colour of the galaxy, and finally, the fit of the 
mass density distribution models (SIS and NFW profiles). 
We have tested it succesfully on simulated data and then, we 
have applied it to a sample of low X-ray luminosity clusters. 

From this analysis we could estimate the mass of seven 
low X-ray luminosity galaxy clusters. One of these clusters 
([VMF981097) was previous ly analysed with a similar ap¬ 
proach (ICarrasco et alj2007f l . We improved the shear fit and 
we obtained a mass consistent with the previous result. For 
the other clusters in the sample, we estimated the mass for 
the first time. 

The velocity dispersions obtained from the SIS fit, are 
in general agreement with the spectroscopic values available 
for four of the clusters in the sample. Masses obtained were 
compared to the X-ray luminosities. Our results are mostly 
in good agreement with pr evious analysis of the M — Lx 
relation, in particular with ILeauthaud et all (l20ld f result. 
In this work we provide further constraints for the M — 
Lx relation, in low-intermediate X-ray luminosity galaxy 
clusters, by increassing the number of observables. 

We plan in future works to include different models to fit 
the shear profiles, in order to include non-spherical models. 
Also, we plan to extend the pipeline to analyse low-massive 
galaxy systems employing stacking techniques. 
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